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We experimentally investigate and analyze the rich dynamics in F—2 spinor Bose-Einstein con- 
densates of 87 Rb. An interplay between mean-field driven spin dynamics and hyperfine-changing 
losses in addition to interactions with the thermal component is observed. In particular we measure 
conversion rates in the range of 10~ 12 cm 3 s~ 1 for spin changing collisions within the F=2 manifold 
and spin-dependent loss rates in the range of 10 -13 cm 3 s _1 for hyperfine-changing collisions. From 
our data we observe a polar behavior in the F=2 ground state of 87 Rb, while we measure the F=l 
ground state to be ferromagnetic. Furthermore we see a magnetization for condensates prepared 
with non-zero total spin. 

PACS numbers: 03.75.Mn, 34.50.Pi, 03.75.Hh 



The investigation of atomic spin systems is central 
for the understanding of magnetism and a highly active 
area of research e.g. with respect to magnetic nanosys- 
tems, spintronics and magnetic interactions in high T c 
superconductors. In addition entangled spin systems 
in atomic quantum gases show intriguing prospects for 
quantum optics and quantum computation Ld- Li LjJ- LJ| ■ 
Bose-Einstein condensates (BEC) of ultra-cold atoms of- 
fer new regimes for studies of collective spin phenom- 
ena [1 0, IIH OH [H HEi BECs with spin degree of 
freedom are special in the sense that their order param- 
eter is a vector in contrast to the "common" BEC where 
it is a scalar. Recent extensive studies have been made in 
optically trapped 23 Na in the F=l state [Tol UlllllllH . 
In addition evidence of spin dynamics was demonstrated 
in optically trapped 87 Rb in the F=l state 0]. There is 
current interest in extending the systems under investi- 
gation to F=2 spinor condensates [ll HE HI HI ITU l20| . 
which add significant new physics. F=2 spinor conden- 
sates offer richer dynamics, an additional magnetic phase, 
the so-called cyclic phase [lJi HH ; as wen as intrinsic con- 
nections to d-wave superconductors |2lj . 

In this letter we present first studies of optically 
trapped 87 Rb F=2 spinor condensates. We measure rates 
for spin changing collisions for different channels within 
the F=2 manifold and discuss the steady state for various 
initial conditions. Additionally we observe and discuss 
the thermalization of dynamically populated mp conden- 
sates. We also present measurements of spin-dependent 
hyperfine decay rates of the F=2 state in 87 Rb, as a key to 
further understanding the intensively studied collisional 
properties of 87 Rb [HE^. 

Our experimental setup consists of a compact double 
MOT apparatus which produces magnetically trapped 
87 Rb Bose-Einstein condensates containing 10 6 atoms in 
the F=2, mp — 2 state. To confine the atoms indepen- 
dently of their spin state they are subsequently trans- 
ferred into a far detuned optical dipole trap. It is oper- 
ated at 1064 nm generating trapping frequencies of typ- 



ically 2-7T x 891 Hz vertically, 27r x 155 Hz horizontally 
and 27T x 21.1 Hz along the beam direction. After trans- 
fer we further cool the ensemble for 500 ms by selective 
parametric excitation resulting in approximately 10 5 
optically trapped atoms and a condensate fraction well 
above 60%. We are able to prepare arbitrary spin com- 
positions using rapid adiabatic passage and controlled 
Landau-Zener crossing techniques [25| at an offset field 
around 25 G. After initial state preparation the magnetic 
field is lowered to a value of typically 340(±20) mG, with 
a field gradient below 15 mG/cm, to ensure a well-defined 
quantization axis and a good overlap of the different mp 
states during spin dynamics |2q|. Spin dynamics is sub- 
sequently allowed to proceed during a variable hold time: 
Stored in the optical trap the condensate spin degree of 
freedom can evolve under well-controlled conditions. Due 
to interatomic interactions, initially prepared mp com- 
ponents can evolve into other mp components, e.g. by 
processes like |0) + |0) «-» |+1) + |— 1), while - disregarding 
atomic losses - the total spin is conserved. This process 
is expected to end in the spinor ground state distribution 
showing the hyperfine dependent magnetic properties of 
the atomic species under investigation. Experimentally 
we detect different spin components by spatially separat- 
ing them with a Stern-Gerlach method during time of 
flight (after switching off the trapping potential). Ab- 
sorption imaging is then used to evaluate the respective 
spatial density distributions as well as the number of 
atoms in the condensate and thermal fractions for each 
spin component. Figure ^ shows typical spinor conden- 
sate evolutions for different starting conditions. These 
pictures demonstrate a rich spectrum of spin dynamics 
in F=2 systems, caused by an intriguing interplay of spin 
component coupling and spin-dependent losses. In the 
following we will investigate these processes separately. 

We analyze the observed s pino r evolution following a 
mean- field approach [lfil IliL Et| |. in which the proper- 
ties of a spinor condensate are determined by a spin- 
dependent energy functional. Extending the approach of 
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FIG. 1: Time-dependent observation of different mp compo- 
nents starting from the initially prepared states denoted by 
(a)-(d). Shown are spinor condensates separated by a Stern- 
Gerlach method (time of flight 31ms). 

|l8| by the experimentally relevant quadratic Zeeman en- 
ergy term we determine the following functional for F=2 
systems: 

K spm = Cl (F) 2 + ^c 2 \(s_)\ 2 - p(F z ) - q(F 2 ), (1) 

where (F), (F z ) and (s_) denote the expectation values 
for the spin vector, its z component and the spin-singlet 
pair amplitude. The experimentally controllable param- 
eters p and q span the phase space for the system. Here 
p represents the mean spin of the system while q is the 
quadratic Zeeman energy, representing the magnetic off- 
set field. 

The spin-dependent mean-field is characterized by the 
parameters ci = ^-^^ and c 2 = 7a "~ 10 7 a2+3a4 , 
introducing the s-wave scattering length a/ for colli- 
sions involving a total spin / of the colliding pair. In 
the following the F — 2,mp = X state is denoted by 
\X). For simplicity relative phases in mixtures are ne- 
glected. The ci term includes all couplings of states 
with Am F = ±1, e.g. |0) + |0) <-> | + 1) + | - 1). 
The C2 term includes the only possible coupling with 
Am F = ±2 : |0) + |0) <-> | + 2) + |-2). The magni- 
tude of these terms is connected to the timescales of spin 
dynamics and their relative strengths indicate the ini- 
tially dominant channels. Taking the calculated scatter- 
ing lengths and our initial experimental conditions (offset 
field B — 340 mG, mean density n « 1.1 x 10 14 cm -3 ) the 
energy ranges for the different terms are: spin indepen- 
dent mean-field ks X 64 nK, spin-dependent mean- field 
ci : fcs x 0..12nK, c 2 : ks x 0..0.05nK and quadratic 
Zeeman energy q : ks X 0..1.6nK. 

Minimizing the energy functional, Eq. leads to the 
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FIG. 2: Spin dynamics and spin-dependent losses, (a) Rel- 
ative population of the mp states of a condensate initially 
prepared in |0). (b) Thermalization of a BEC in the |+2) 
state form initially |0) + |+1). (c) Total spin of the condensate 
initially prepared in | + (d) Decay curves for total number 
of condensed atoms in all mp states, given for different inital 
m,F preparations. 

parameter-dependent ground state spin composition of 
the system, reflecting its magnetic properties 28] . The- 
ory predicts 87 Rb in F=2 to be in the polar phase, but 
close to the border to the so-called cyclic phase [l^ . 

Experimentally we find the system dominated by three 
processes, setting different timescales, in increasing or- 
der: spin dynamics, two-body hyperfine losses and three- 
body recombination. This hierarchy allows us to sepa- 
rately analyze individual interaction processes, which are 
exemplarily given in Fig. [21 

In a first set of measurements we study the evolution 
of spinor condensates starting from different initial spin 
states. Figure^a) shows the evolution of the relative oc- 
cupation of the different mp components of a F=2 spinor 
condensate initially prepared in the state |0), also shown 
in Fig. Ufa). This example clearly demonstrates that spin 
evolution takes place predominantly between neighboring 
spin states, as expected from coupling via the dominat- 
ing ci-term discussed above. Thus first the components 
|±1) and only afterwards the |±2) states are populated. 
Interestingly there is a delay in the occupation of |±1) 
and a further delay in the occupation of |±2) states. The 
latter delay can be intuitively explained since there must 
first be some occupation in the |±1) states before the fully 
stretched states can be occupied. The origin of the initial 
delay as well as the "overshoot" of the |±1) occupation 
is still under investigation. 

As an interesting feature Fig. E^a) shows evidence of 
spinor oscillations on a 100 ms timescale around equipar- 
tition. This indicates a coherent evolution between the 
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TABLE I: Spin evolution for different preparations. 



initially 


initial 


initial channels 


Anally 


prepared 


total 


into m,F state 


populated 


niF states 


spin 


G [lO-^crn^- 1 ] 


m,F states 


10} 





-» |±1) « 21.0 


equipartition 


| + 1) + |-1) 





-> |0) » 26.9 


equipartition 






-> |±2) « 4.6 




1+1) + 10) + 1-1) 





-> |±2) w 5.0 


equipartition 


|+2) + |-2) 







1 + 2) + 1-2) 


| + 2) + |0) + |-2) 





-!■ |±1) < 0.1 


1 + 2) + 1-2) 


1 + 2) + 1-1) 


1/2 




1 + 2) 


1 + 1) + 10) 


1/2 


| + 2) « 21.7 


1 + 2) 






-> |-1) « 19.2 




1 + 1) 


1 


-> j + 2) « 22.4 


1 + 2) 






-> |0) » 12.2 








H|-l)«4.7) 




1 + 2) 


2 




1 + 2) 



different mp components, with only minor influence of 
the thermal cloud on this timescale. 

In fact we always observe that spin dynamics is sig- 
nificantly quicker than thermalization. Spin dynamics 
thus leads to the formation of pure condensate wavefunc- 
tions in new spin states, which subsequently thcrmalizc 
on a much longer timescale. This opens a new way to 
study thermalization effects with an independently tun- 
able heat bath. In Fig.^b) the condensate and thermal 
fractions of the newly created |+2) state, starting with a 
mixture of the states | + 1) and jO) are plotted. At first 
the condensate in the | + 2) state is almost pure. This 
can be intuitively understood, as spin dynamics predom- 
inantly occurs in the high density condensate fraction 
on a timescale of the order of 25 ms. The subsequent 
buildup of a thermal cloud in the newly populated spin 
states takes place within a significantly longer time, of 
the order of 100 ms. Dominant effects contributing to 
the thermalization of the new spin state are the interac- 
tion of the new spin condensate with the parent thermal 
cloud ("melting"), and spin dynamics of the parent con- 
densate fraction with the parent thermal cloud leading to 
a direct occupation of the new component thermal cloud. 
In this respect spinor condensates are a possible new tool 
for investigations of finite temperature effects. 

In order to draw a more complete picture of spin dy- 
namics in F=2 systems we now investigate the evolution 
of a variety of initial states as summarized in Table [I] In 
this table we give estimates for the initial spin dynamics 
rates for different channels as well as the experimentally 
observed "final" distribution. The rates are intended for 
comparison with loss rates and between the different spin 
channel rates. They were obtained using the total ensem- 
ble density, number of BEC atoms in the prepared mj? 
components and the initial slope of a phenomenological 
exponential fit to the respective spin component popula- 
tion curves: 



At short timescales we find the rate for |0) — ► |±1) similar 
to the rate |±1) — ► |0) as one would expect for a reversible 
or coherent process. In contrast, at longer timescales the 
system evolves as expected into a time-independent dis- 
tribution of rciF states. As a main result our measure- 
ments show the stability of a mixture of |±2) (see Tabled 
and Fig. [Jd)), for which we also find no spatial separa- 
tion of the two components in the trapping volume. This 
is clear evidence of polar behavior as predicted for F=2 
spinor condensates of 87 Rb. 

Some initial preparations with zero total spin end up 
with all rriF components equally populated (see also 
Fig. ^a,c)). This is in fact one of several degenerate 
ground states of the polar phase in the case of zero mag- 
netic field |20j, but the quadratic Zeeman energy lifts 
this degeneracy. For our experimental parameters a mix- 
ture of all to f components is not a ground state for any 
phase in the mean-field description. The observed behav- 
ior might be due to a lack of time to reach the ground 
state, as spin-dependent losses depopulate the conden- 
sate before thermal equilibrium is reached. 

Remarkably, if the atoms are prepared in superpo- 
sitions which are ground states for the cyclic phase 
we observe very slow dynamics. First the mixture of 
|+2) + |0) + | — 2), which is the ground state of the cyclic 
phase with total spin zero |20| shows nearly no spin dy- 
namics. If starting with 1+2) + |— 1) as a ground state for 
the cyclic phase at B > |2£| with non-zero spin we do not 
observe any spin dynamics. These cases are particularly 
surprising, as the states |0) alone as well as | + 1) alone 
show fast spin dynamics as can be seen in Fig.^a,b). In 
the mixtures however, we observe only a faster decay of 
the |0) and |— 1) component compared to |±2). Conclud- 
ing this, the stability of the |±2) mixture demontrates 
polar behavior of F=2 87 Rb atoms, and in addition the 
slow dynamics of prepared cyclic ground states shows the 
F=2 state to be close to the cyclic phase Ejj . 

An important point in the dynamics of spinor con- 
densates is total spin conservation. This is directly ob- 
served in the rates for all spin preparations (see rates 
Table HJ) except in the case starting with | + 1). Here 
we observe a much higher rate for the production of the 
|+2) atoms than for the ones in |0). This can be under- 
stood as follows: Imagine the |0) atoms favor to change 
their spin into — 1) while interacting with an atom in 
the initial spin state which then changes into |+2). Then 
the correct rates are determind by the measured values 
G^|+2) — G->\-i) an( i ^^l-i) • These two values 

are 17.7 x 10 _13 cm 3 s _1 and 16.9 x 10 -13 cm 3 s -1 , respec- 
tively, and agree within the expected precision. 

Investigating spinor condensates with initial spin ^ 
one finds that the combination of spin dynamics and spin- 
dependent hyperfine losses (see below) leads to a loss 
induced magnetization during the evolution to the final 
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state (Fig. njb)[5Jc)). The remaining condensate always 
ends up in the fully stretched state if the initial state was 
prepared with non-zero total spin or in any other non- 
symmetric mixture. This is due to the fact that only 
the fully stretched component is immune to hypcrfmc- 
changing collisions. The magnetization process is only 
inhibited for a symmetric initial state having total spin 
zero. 

Next we focus on hyperfine-changing collisions (i.e. in- 
volving transitions F=2— >F=1) as a special case of spin 
relaxation dynamics. The release of hypcrfine energy in 
these collisions leads to immediate loss of the collision 
partners from the trap. This loss process usually takes 
place on very short timescales, prohibiting the observa- 
tion of spin dynamics in the upper hyperfine level, e.g. of 
23 Na [15j. The relatively low hyperfine loss rates we find 
for 87 Rb are due to a coincidental destructive interference 
of decay paths 0,H3|- Fig.EJd) shows the decay of Bose- 
Einstein condensates of 87 Rb prepared in different mp 
states/mixtures. The decay of the stretched state |+2) 
is dominated by three-body recombination and was mea- 
sured by to be L = 1.8 x 10~ 29 cm 6 s _1 , a value which 
also fits to our data. The equal superposition of the |+2) 
and I — 2) states is special in the sense that we observe 
no spin dynamics in this case. The hyperfine-changing 
collisions can thus only occur via collisions of the type 
|2,+2) + |2,-2) -> |l,mi) + \F,m 2 )(F = 1,2). We de- 
duce the two-body rate G = 6.6(±0.9) x 10 _14 cm 3 s _1 , 
which leads to a much faster decay than three-body losses 
for the pure |+2) state (Fig.Efd)). The dependence of the 
hyperfine-changing collision rate on the number of decay 
channels involved becomes obvious observing the decay 
of a condensate initially prepared in the state |0) (see 
Fig-EJd)). It is important to note that in this case spin 
dynamics more rapidly (< 25 ms) distributes the popula- 
tion almost equally over all nip states. The decay curve 
(of initially |0)) thus effectively represents the loss of a 
spin state equipartition, which has access to all possible 
channels. We determine the mean two-body decay rate 
for this case to be G = 10.2(±1.3) x lO^cmV" 1 . 

Finally we have also studied spinor condensates in the 
F=l state where we observe slow spin dynamics on a 
timescale of seconds like in the case of 23 Na in F=l [Icj . 
Starting in the superposition state |1,±1) we observe 
the creation of the |1,0) component. The final state 
with all mp components populated and spacially mixed 
is reached after about seven seconds. According to the 
phase diagrams of we have thus shown that 87 Rb in 
the F—l state is ferromagnetic as predicted by 

In conclusion we have presented studies of various as- 
pects of F=2 spinor condensates in 87 Rb. In particular 
we have discussed rates for spin dynamics as well as the 
evolution into the magnetic ground state. Furthermore 
we have investigated spin-dependent hyperfine-changing 
collisions and thermalization effects in newly created spin 
components during spinor evolution. These studies have 



demonstrated the feasibility of 87 Rb condensates as a 
model for F=2 spin systems, added data to the collisional 
properties of 87 Rb and opened a new path for finite tem- 
perature studies. As a key result we measure 87 Rb in the 
F=l state to be ferromagnetic and observe polar behav- 
ior for the F=2 state. 

We acknowledge stimulating contributions by W. Ert- 
mer and support from the Deutsche Forschungsgemein- 
schaft in the SFB 407 and the SPP 1116. 
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